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THE ANNULAR MOMENTUM CONTROL DEVICE (AMCD) 

AND POTENTIAL APPLICATIONS 

By Willard W, Anderson and Nelson J. Groom 
Langley Research Center 

SUMMARY 

An annular momentum control device (AMCD) consisting principally of a spinning 
rim, a set of noncontacting magnetic bearings for supporting the rim, a noncontacting 
electric motor for driving the rim, and, for some applications, one or more gimbals is 
described. The device is intended for applications where requirements for control torque 
and momentum storage exist. Hardware requirements and potential unit configurations 
are discussed. Theoretical considerations for the "passive” use of the device are dis- 
cussed. Potential applications of the device in other than passive configurations for the 
attitude control, stabilization, and maneuvering of spacecraft are reported. 

INTRODUCTION 

The use of stored angular momentum for the purpose of controlling the attitude of 
fine pointing or long duration spacecraft, where either environment contamination or 
excessive fuel use prohibit reaction jet usage, has become almost universal. Applications 
of the concept include spinning spacecraft (Tiros), dual -spin spacecraft (OSO), momentum 
wheel stabilized spacecraft (ITOS), reaction -wheel stabilized spacecraft (OAO), and, 
recently, the control moment gyroscope system used for stabilizing Skylab. 

Each of these methods has advantages and disadvantages. Spinning the spacecraft 
(or a portion of it) to achieve attitude stability is simple and reliable, yet the spacecraft 
itself cannot be utilized fully because of its rotation. Also, any momentum axis reorien- 
tation maneuvers require external torques for momentum precession and the artificial 
gravity field produced may run counter to payload requirements for zero gravity. Sta- 
bilizing the spacecraft by utilizing a momentum wheel which provides gyroscopic stiffness 
equivalent to spinning the vehicle itself allows a nonspinning spacecraft and permits arbi- 
trary orientations about the roll axis for the purpose of pointing onboard experiments. 

This technique does not overcome the inability to reorient or maneuver the spacecraft 
about all three axes without external torques since the spin axis of the momentum wheel 
is fixed with respect to the spacecraft. The use of three reaction wheels with axes alined 
with the spacecraft axes allows complete spacecraft active attitude control. However, 



reaction -wheel momentum must be limited to relatively low amounts because of an 
excessive requirement for power when directly producing a torque on a rapidly spinning 
flywheel. The limitation on react ion -wheel momentum can be overcome by using a con- 
trol moment gyro (CMG) system, which uses constant speed wheels and develops preces- 
sion torques through controlled slewing of gimbaled flywheels. However, to achieve the 
smooth low-level torques necessary for fine pointing requires precise control of very 
low gimbal slew rates. These low gimbal rates are inherently limited by the requirement 
for high servo stiffness and thus high friction torque, unless extreme mechanical preci- 
sion and resultant high cost are involved. 

The purpose of this paper is to introduce a new development in the field of momen- 
tum storage devices, the annular momentum control device (AMCD).* The spinning 
assembly of the device is discussed in detail and its advantages relative to conventional 
flywheel design are presented. The use of the spin assembly as a momentum wheel 
('’passive" case) with spin axis control for spacecraft with orientation requirements such 
as those of Tiros, ITOS, or TACSAT is described in detail. 

The use of the device for active pointing control, maneuvering, and other applica- 
tions of spacecraft and experiment control is discussed in general terms. Emphasis is 
placed on the use of the device in applications where a large radius is allowable and ben- 
eficial since those applications are new. The device can be used as the spin assembly 
for any of the conventional momentum storage devices including reaction wheels and con- 
trol moment gyros. 


SYMBOLS 

Values are given in both SI and U.S. Customary Units. The measurements and 
calculations were made in U.S. Customary Units. 

A rim cross-sectional area 

E transformation matrix 

F magnetic bearing forces 

G external torques 

H angular momentum 


^Patent pending for "Annular Momentum Control Device," NASA Case 
No. LAR 11051-1. 
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H a 

H s 

I 

la 

Is 

Im( ) 
K £ 
Kfi 
K s 
K5* 

K X 

K X 

k 

m 

Re() 

r 

s 

T 

T R 
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total angular momentum of AMCD 
spacecraft total angular momentum 
moment of inertia 
transverse AMCD rim inertia 
transverse spacecraft inertia 
imaginary part of complex root 
magnetic bearing linear displacement gain 
magnetic bearing linear rate gain 
flywheel energy shape factor (see eq. (Al)) 
flywheel momentum shape factor (see eq. (1)) 
magnetic bearing angular displacement gain 
magnetic bearing angular rate gain 
radius of gyration 
mass 

real part of complex root 
radius 

distance along rim; also, Laplace operator 
interaction torque 
flywheel energy 
time 


3 



V 

X,Y,Z 

Y 
a 
0 
r 
6 

e H 

e I 

n 

9 

X 

« 

P 

(Psyst) max 

Poo 

a 

T 

<P 


arbitrary vector 
coordinate axes 
Young’s modulus 
telescope azimuth 
telescope elevation 

angle between axis of rotation and X-axis (see eqs. (C14)) 
spacecraft rotation angle 
spin momentum ratio, H s /H a 
transverse inertia ratio, I a /I s 

rim vibration amplitude 

transformation angle, Y-axis rotation angle 

rim rotation angle about axis of rotation (see eqs. (C14)) 

rim rotation angle 

mass density 

maximum system damping (see eq. (12)) 

damping coefficient 
working stress 

magnetic bearing gain ratio, K^/K^ 

X-axis rotation 
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uncoupled rim precessional frequency 


« angular velocity 

<*>' angular velocity cross product matrix 

rim fundamental frequency 

Wjj approximate precession frequency, (H a + H s )/(I a + I s ) 

Woo stiffness factor, (K^/I s )V2 

Subscripts: 

a AMCD rim 

b body 

d damped 

f body fixed 

i inertial 

s spacecraft 

T total 

x,y,z coordinates 

Superscripts: 

c coordinate axes 

-1 matrix inversion 

* nondimensional with respect to u>jj 

Dots over symbols denote derivatives with respect to time. 
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AMCD SPIN ASSEMBLY DESCRIPTION 


The spin assemblies of momentum wheels, reaction wheels, and control moment 
gyroscopes used in attitude control systems for spacecraft to date have been of nearly 
conventional terrestrial design (shaft -driven steel flywheels with ball bearings). The 
AMCD spin assembly configuration is based on space usage (vacuum and zero gravity) 
and consists of a rotating rim (no central hub) suspended by noncontacting magnetic bear- 
ings and powered by a noncontacting linear electromagnetic motor. (See fig. 1.) 


Magnetic bearings 



The choice of configuration is based on several factors. The first factor is a 
desire to obtain maximum angular momentum for a given flywheel mass. The relation- 
ship of these two variables is discussed in appendix A (eq. (A3)) and is given as follows: 



where H is the angular momentum, m the mass, k the radius of gyration, K s a 
dimensionless shape factor, cr the allowable working stress, and p the density of the 
flywheel. A consideration when optimizing H/m is the specific stress a/p of the 
material. From reference 1, working values of the square root of specific stress for 
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conventional metals (steels and aluminums) lie below 350 m/sec (1150 ft/sec), the 
stronger maraging steels approaching 500 m/sec (1640 ft/sec). Values of specific stress 
for the filaments of the stronger composite materials range from 1060 m/sec (3480 ft/sec) 
for boron and graphite to 1230 m/sec (4050 ft/sec) for the Dupont material PRD-49 
(Kevlar). Clearly, from the aspect of specific stress, the composite materials are better. 
However, to achieve full utilization of the material strength requires a unidirectional 
layup of the composite. 

Another consideration is the combined effect of size (radius of gyration) and shape. 
For a flywheel of constant shape, H/m is proportional to the radius of gyration k, and 
it is clear, therefore, that the largest possible k is optimum. However, flywheel size 
is normally limited by outside diameter and not by the radius of gyration. Equation (A3) 
is thus rewritten as 



where r Q is the outside radius of the flywheel and 



The maximum value for Kg for the nine flywheel shapes considered in reference 2 
is Kg = 0.500 for the thin rim. It is interesting to note that not only is the thin rim an 
optimum flywheel shape but it is also one of the few flywheel shapes that will allow full 
utilization of the filament strengths of composite materials by allowing a unidirectional 
layup. This fact also makes the AMCD an attractive energy storage device. 

The second factor behind the selection of the AMCD configuration is the multiple 
problems of wear, reliability, and vacuum lubrication of the ball bearings used in current 
and planned control moment gyros. To achieve the large momentum needed for Skylab 
and for planned shuttle and shuttle -launched missions requires flywheel speeds at values 
where ball bearing design is an empirical art necessitating extensive preflight testing of 
the bearings and lubrication system. As has been evidenced on Skylab, where two out of 
three of the control moment gyros developed bearing problems (one causing complete 
failure of the unit), poor performance can still result from this empirical design process. 
The use of magnetic bearings to support the rim eliminates all wear, requires no lubri- 
cation, and increases the reliability to that of the solid-state electronic bearing circuits. 

A third factor is the isolation of the rotating rim from the spacecraft afforded by 
the magnetic bearings. Spacecraft vibration resulting from reaction wheel or CMG rotor 
unbalance and ball bearing element rotation can have a detrimental effect on the fine 
pointing of vehicles with high accuracy requirements such as the planned large space 
telescope. (See ref. 3.) Although imperfect balance of an AMCD rim would cause rim 
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motion, for the case of all active nonpermanent magnet bearings the motion itself would 
produce no disturbance on the spacecraft unless detected by the bearing sensors and then 
acted upon by the bearing force coils. For this reason, appropriate notch filtering at the 
spin frequency and at the fundamental vibrational frequency and possibly at the higher 
harmonics of both should eliminate vibrational input to the vehicle. It is interesting to 
note that for the thin -rim AMCD, it is not possible to have resonance between spin and 
spin plane vibration frequencies because the rim vibration frequencies, for the assump- 
tions of appendix B, are always higher by a factor of at least two. Although the assump- 
tions of appendix B do not cover all potential AMCD rim geometries, the equations illus- 
trate the stiffening which occurs because of the tangential stress generated by rim 
rotation. Also this stiffening has the effect of circularizing the rim and improving bal- 
ance at design speed. 

A final factor to consider is the ability to command output torques for control 
usage directly with no mechanical or electrical breakout torque. This ability is implicit 
when the device is used either in the ’’passive’' mode to be discussed later where control 
torques are generated by a damped centering of the rim between bearing elements or by 
direct command from spacecraft sensors when used in the "active’’ mode which is also 
discussed later. 


AMCD SPIN ASSEMBLY HARDWARE CONSIDERATIONS 

Because of its unique configuration, the AMCD presents several design constraints 
which interact with all components and result in basic requirements for any AMCD 
design. Obviously, any number of variations in the design of rim, motor, and/or bearings 
could be made in an attempt to better meet the specific design requirements for a given 
mission, but some general considerations warrant discussion. 

Magnetic Bearing Considerations 

For AMCD designs of large radial dimension, a very large weight penalty could 
occur if a continuous magnetic bearing around the rim was provided. In practical terms 
the total magnetic bearing structure weight could be two to three times that of the rim. 
This condition could compromise any system weight advantage presented by the AMCD 
design. The practical solution to the bearing problem is to use segmented bearings as 
shown in figure 1. The minimum number of bearing segments would be three with a 
variable number of elements per segment depending on particular design requirements. 

Unfortunately, the use of segmented bearings to solve the bearing weight problem 
produces another potentially serious problem. Since the bearings are now in discrete 
sections, the rotating rim will see a discontinuous magnetic field. Unless proper choice 
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of rim design is made, this could result in appreciable eddy current and hysteresis 
power losses in the rim portion of the bearing magnetic circuit. 

Rim Considerations 

As mentioned previously, the use of segmented magnetic bearings introduces the 
design requirement that the rim portion of the bearing magnetic circuit have low eddy 
current and hysteresis losses. Thus, a low reluctance material with a very high resis- 
tivity and a narrow hysteresis loop is required. A material that has these characteris- 
tics in addition to the structural characteristics required for the rim structure does not 
currently exist. A practical solution to the conflicting rim material problem is to use a 
high-strength nonmagnetic and nonconducting composite material such as graphite epoxy 
for the rim structure and embed a low-loss magnetic material such as a ferrite material 
in the main rim structure to serve as the rim portion of the bearing magnetic circuit. 

Rim Drive Motor Considerations 

As with the magnetic bearings, the weight penalty for a motor stator that is contin- 
uous around the rim would be too high. Use of a linear motor employing one or more 
stator segments is one obvious solution. Again the problem of conflicting requirements 
for the rim material arises. As in the magnetic bearing solution, the rotor part of the 
motor magnetic circuit can be embedded in the main rim structure. For example, a 
permanent magnet dc brushless drive motor could be designed by utilizing permanent 
magnets embedded in the rim with a magnetic field detector to sense the position of a 
given magnet for commutation purposes. 

AMCD Laboratory Model 

An AMCD laboratory model embodying the design philosophy presented is being 
built for the Langley Research Center by Ball Brothers Research Corporation. The rim 
(fig. 2) is made of a graphite-epoxy composite, is 1.7 meters (5.5 feet) in diameter, 
weighs 222.4 N (50 pounds), and will rotate at 3000 rpm. The suspension system con- 
sists of three segments spaced equidistantly around the rim. Magnetic bearing elements 
at each segment interact with low-loss ferrite material embedded in the rim. The rim 
will be spun up by a series of stationary electromagnets that both push and pull on small 
samarium-cobalt permanent magnets also embedded in the rim. 

THEORETICAL CONSIDERATIONS OF A "PASSIVE” 

AMCD -SPACECRAFT SYSTEM 

The use of the AMCD as a "passive" device involves rigidly mounting the magnetic 
bearings and spin motor to the spacecraft and providing an active spacecraft control loop 
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Figure 2.- Graphite -epoxy rim for the AMCD 
laboratory model. 


using spin motor torques for spacecraft roll or spin control. (See fig. 3.) The word 
passive is used because spacecraft sensors are not required for stabilization of space- 
craft pitch and yaw axes although rim position sensors are required for the active mag- 
netic bearings. Nonlinear and linearized equations of motion for the passive AMCD- 
spacecraft case are derived in appendix C. These equations were used to generate time 
histories and root locus plots which are discussed later. Before presenting these data, 
the rotational stability of the device is discussed. 

Basic Rotational Stability Considerations 

The linear model derived in appendix C does not contain spin axis torques when the 
AMCD spin motor torque is controlled to counteract AMCD rim friction (hysteresis and 
eddy current losses) since both the AMCD rim and the spacecraft are assumed to be sym- 
metrical about the spin axis. This means that there is no mechanism within the model to 
allow the coning angle to increase since any increase in the coning angle can only be the 
result of a decrease in spin momentum if there are no external torques present. 
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Figure 3.- ’’Passive" AMCD-spacecraft system. 


To explain this further, assume that the relative transverse rotation between rim 
and spacecraft is small (as has been assumed throughout the report). Thus 

Ht cos 0 = ^az^az + ^sz^sz (2) 

where Hq< is the total angular momentum of the system, 9 is the angle between the 
Z -axis and the momentum vector, I az and I sz are the rim and spacecraft spin-axis 
moments of inertia, and cu az and o> sz are the spin rates of the rim and spacecraft. 
Differentiating with respect to time yields 

a _ laz^az +Isz^sz 

H t sin 9 W 

and thus an increase in the cone angle 6 is seen to require a net decrease in spin 
momentum. Obviously, if there are no spin axis torques acting directly on the two bodies 
and the two bodies are symmetrical about the spin axis (no spin -axis cross-coupling 
torques), the spin momentum is constant and therefore 9 is constant. 

In any real physical system, internal sources of energy dissipation do exist which 
can change the spin momentum. There are many references in the literature dealing with 
the subject of rotational stability based on models of dual-spin vehicles containing internal 
sources of energy dissipation which affect spin momentum. Vernon Landon (ref. 4) is 
credited with being the first investigator to explain the stability of dual -spin vehicles when 
these vehicles are not spinning about an axis of maximum inertia. For the purposes of 
this paper, it is assumed that the energy dissipation sources exist almost entirely within 
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the vehicle rather than within the AMCD rim. For this reason the requirement for rota- 
tional stability (for the case of small relative transverse rotations) reduces to a require- 
ment that the rate of precession of the spin axis (Z-axis) about the total momentum vector 
be greater than the rate of spin of the spacecraft. (See ref. 5 for a mathematical deriva- 
tion of this requirement.) For the passive AMCD-spacecraft system, this condition 
represents a minimum requirement for AMCD rim momentum which must be met to 
insure asymptotic rotational stabilities. 

The linear model does contain terms which can render the system unstable. How- 
ever, the instability would not affect the system cone angle, but rather would allow the 
relative transverse rotation between rim and spacecraft to increase in an unstable man- 
ner. These instabilities are examined after the characteristic equation for the linear 
model is discussed. 

Characteristic Equation for the Linearized Model 

The nonlinear and linearized equations of motion for the passive AMCD-spacecraft 
case are derived in appendix C. The nonlinear equations were used to generate time 
histories which are discussed later. The linearized set of equations (eqs. (C21)) has 
the following characteristic equation which will be used to discuss system behavior: 
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and where and are the resultant angular and angular rate magnetic bearing 

constants, I a and I s are the AMCD rim and spacecraft transverse moments of inertia, 
and H a and H s are the AMCD rim and spacecraft spin momentums. (See appendix C.) 

The coefficients of equation (4) were used to form a symbolic Routh stability array 
by means of a computerized symbolic manipulation routine called MACSYMA. An attempt 
to derive simplified conditions for system stability was not successful because of algebraic 
complexity. The stability of the system was then examined by using the method of refer- 
ence 6. This examination yielded the following two conditions for stability for positive 
damping r: 


1 + e H 



1 + Cj 

i 1 + e nf 




c 

e H R 



>0 


These conditions agree numerically with the Routh array conditions. 

It is informative to discuss simplified versions of equation (4) initially. For the 
case where the AMCD rim transverse inertia I a and the spacecraft spin momentum H s 
are relatively small (ej and ejj — 0), equation (4) simplifies to 


|( T 2 


Woo 4 + 


W H 2 )s 4 + (2TW.X, 4 + 2tw 0O 2 w h 2 )s 3 + (wco 4 + T 2 U) co 4 W H 2 + 2w 00 2 w h 2 )s 2 


( 2ra, oo 4a, H^) s + (* 


W 


H 


= 0 


(6) 


Referring to the first of equations (5), the variable r is seen as a ratio of mag- 
netic bearing gains or the required lead on the magnetic bearing rim displacement sig- 
nal. (Also see eqs. (C17).) This variable defines system damping. Insight into the 
variable can be gained by letting the AMCD momentum become large {w R - <*>). 

Equation (6) is reduced to 

S 2 (s 2 + TWo^S + Wo^) =0 (7) 

The combined behavior described by this equation is that of uncoupled X- and Y-axis 
vibration of the spacecraft (the two complex conjugate roots) with damped frequency 



and damping coefficient 
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and precession of the system (the two s = 0 roots) about the X- and Y-axes caused by 
external torques about the Y- and X-axis, respectively. Since the momentum is large 
for this case, the rate of precession (which is inversely proportional to momentum) 
would be zero in the limit as approached infinity. 

The variable Wjj can be defined by allowing the control gain ratio or lead t to 
approach infinity. Under this circumstance the AMCD rim is not allowed transverse 
motion with respect to the spacecraft; this situation is thus likened to a nonflexible dual- 
spin spacecraft. Equation (4) reduces to 


s 2 ( s 2 + 5 2 ) (s2 + w H 2 ) = 0 (8) 

Thus, the variable Wjj is the frequency of the oscillatory mode of the system when no 
lelative transverse motion is allowed. It is noted that the system has no damping when 
r approaches infinity. Since the system obviously has no damping when r is zero, the 
selection of the value of r to maximize system damping is critical. 

Figure 4 contains root locus plots for equation (6) with respect to a set of nondi- 
mensional variables, 


* - s _ 


s* - 


CO 


H 


= Re* + i Im* 




T* = TC 0 


H 


* - Woo 


WJ, = 


w 


H 


(9) 


This relation allows a single figure representation, and also from practical considerations 
w R is of the order one. Only the upper left quadrant of the s* plane is presented, the 
lower quadrant being a mirror image. Values of to* = 1.0, 1.5, and 2.0 and t* vary- 
ing from zero to infinity were used to generate figure 4. 


For no damping (t* = 0), the two positive roots of equation (6) are 


s i,2 = ia, » 




+ 1 ± w* 


(w * 2 + 4) 


1/2 


1/2 


( 10 ) 


If the damping t* is increased, the root loci are seen to be approximately semi- 
circular, the positive limiting roots defined by equation (8) for r* large being 

s* = 0, 0, i (n) 

This situation corresponds very nearly to that for a spin -stabilized, dual-spin- 
stabilized, or momentum-wheel-stabilized spacecraft and points out a major advantage of 
an AMCD passively” stabilized spacecraft, namely, the introduction of system damping. 
The maximum value for this damping for the ” system” mode (the root which approaches 
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s * - i or s = i^jj) which defines the "coning" motion associated with spin -stabilized 
spacecraft can be approximated by assuming the root locus to be semicircular. This 
assumption will yield 


( p syst) max 



1 + s* 2 


( 12 ) 


Although equation (12) predicts lower damping than that seen in figure 4, it is sufficient 
to illustrate trends. Figure 5 is a plot of (p s yst) max determined from equation (12) 
for selected values of co,*, and u> H . As can be seen from the figure, increasing the 
AMCD momentum (increasing u> H ) or decreasing the bearing stiffness (decreasing a> M ) 
will allow higher maximum system damping. 



From figure 5 it is seen that values of (Psyst) max could theoretically be as high 
as one. However, from hardware considerations, there is a practical limit to the lead 
r which must be taken into account in any specific design. An optimum selection of 
system damping has been accomplished which minimizes an integral absolute spacecraft 
pointing error criteria following an initial spacecraft offset. The value for the dimen- 
sionless variable t* based on these simulation runs is 


r* = 2.8 
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for a selected value of of 1.0. This value of r* was independent of o>jj for the 

values of Wpj selected which varied from Wjj = 0.25 to Wjj = 1.0. This value is 
circled in figure 4 and can be used as a basis for an initial design of an AMCD ’’passive” 
system, provided ej and ejj are negligible. 


Spacecraft Rotation Effects 

For some applications nonnegligible spacecraft rotation may be a requirement 
( e H * 0) and would occur, for example, for missions where Earth or satellite tracking was 
necessary. Figure 6 contains root locus plots for increasing r* for values of 6 jj = 0, 
0.25, and 0.5; = 1.0; R = 2; and ej = 0. Since the variables plotted are normalized 

with respect to u> H , the plots can be considered to represent a system with fixed total 
momentum. Therefore, for a given vehicle, transferring momentum from the AMCD to 
the spacecraft for the case where the spacecraft and AMCD rotate in similar directions 
( e H > o) can be represented by loci of roots of constant t* and increasing ejj. The 
’’system" root exhibits decreased damping resulting from the decrease in AMCD momen- 
tum and this effect would require modification to these design criteria (r* = 2.8). 

Figure 7 contains root locus plots for increasing t* for values of ejj = 0, —0.1 , 
and -0.15; 1; R = 3; and €j = 0. Again, a given vehicle can be represented by loci 

of roots of fixed r* as ejj increases negatively. The system root exhibits increased 
damping, resulting from increased AMCD momentum. 


AMCD Inertia Effects 

The AMCD rim has a basic precessional mode with an uncoupled frequency (for the 
thin rim assumption, I az = 2I a ) of 


H 


«p = ~ = U 

■la 


Since the rim rotation rate | is large relative to other system frequencies, the 
effect of the rim precessional mode on the system roots has been found to be negligible; 
for example, in the generation of figures 4 to 7, no noticeable effect occurs when antici- 
pated values of €j are included, and roots are found by use of equation (4). However, an 
additional complex root is found which represents rim precession, and which requires ade- 
quate magnetic bearing damping (proper selection of r*), Figure 8 contains plots of the 
damping ratio for this rim mode as a function of bearing lead r* for values of 6pj = -0.5 
to €jj = 1.0 and for a value of 6j = 0,0025 which results when nominal parameters are 
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Figure 8.- Rim damping ratios. 6j = 0.0025. 

selected. Values of the damped natural frequency of this mode did not vary significantly 
from the uncoupled precessional frequency J2p. 

Simulation Results and Discussion 

The nonlinear equations of motion for a passive AMCD-spacecraft configuration, 
which are derived in appendix C, were used to provide the time histories presented in 
this section. The spacecraft that was simulated does not represent a particular vehicle 
but should be representative of a typical Earth observation satellite. The configuration 
is basically a right circular cylinder with a diameter of 1.52 meters (5 feet), a length of 
2.5 meters (8.2 feet), and a mass of 1175 kilograms (80.5 slugs). The transverse iner- 
tias are 680 kg-m^ (500 slug-ft^). An experiment-induced disturbance was used to 
illustrate the short-term dynamics of the system. The disturbances, shown in figure 9, 
represent the operation of protective doors and were taken from reference 7 (p. DC-17). 
The AMCD parameters were derived by assuming a thin -rim configuration with a diameter 
of 1.52 meters (5 feet) and a tip speed of 305 m/sec (1000 ft/sec). The parameters for 
the magnetic bearings (that is, stiffness and rate gain) were derived from the information 
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presented in figure 4 for o> H and of 1. From equations (5), since and e H 

approach 0, the momentum H a of the AMCD is equal to the transverse inertia I s of 
the spacecraft. From equations (9) (with cujj = l), t* = r and cu& = w*,. By using the 
expression for oj^ 2 in equations (5), the magnetic bearing stiffness was found to be 
7.78 N/cm (4.44 lb/in.). The r* for optimum system damping, as indicated in figure 4, 
is 2.8, For purposes of illustration, one value of t* above and one value below the r* 
for optimum damping were included in the simulation. These values were 6.4 and 0.8, 
respectively. 



Figure 9.- Experiment disturbance torques. 


The results of the simulation runs are presented in figures 10 to 14. Figures 10 
to 12 represent a system with r* = 2.8. Figures 10 and 11 show system pointing accu- 
racy and damping characteristics. The system damping is very good, system motion 
being essentially damped out in approximately 3 cycles. Figure 12 is presented for the 
purpose of illustrating practical air gap requirements. It represents the motion of the 
rim center point within particular axial magnetic bearing segments so that the total peak- 
to-peak motion would represent the total bearing gap required to allow uninterrupted 
motion of the rim. The minimum bearing gap for the case considered is approximately 
0.0012 cm (0.00047 in.). The bearing gaps would be made much larger than this in prac- 
tice in order to ease machining tolerances. Figures 13 and 14 show system pointing 
accuracy and damping characteristics for r* of 0.8 and 6.4, respectively. 

DESCRIPTION OF OTHER POTENTIAL USES OF AN AMCD 
An "Active” AMCD -Spacecraft System 

The use of an AMCD as an active device involves the generation of spacecraft 
torques in direct response to sensed spacecraft pointing errors. The torques would be 
generated by energizing the AMCD magnetic bearing flux coils in response to desired 
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Time, sec 

Figure 14.- Response of spacecraft to disturbance input with 
u>£, = 1 and t* = 6.4. 

control torques based on spacecraft sensor outputs. The spacecraft is reoriented by 
producing torques against the AMCD rim momentum. 

Two possible configurations of an active system are shown in figure 15. The spin- 
axis -pointing configuration of figure 15(a) allows both spacecraft pitch and yaw axes to 
be controlled by the magnetic bearing/torquers. However, if the target were the Earth 
or Sun, this configuration would require application of an external torque to precess the 
AMCD momentum vector at the appropriate rate. This configuration might be practical 
for solar viewing from Earth orbit (1° per day precession) by producing torques against 
the Earth's magnetic field with the use of conventional flux coils. Figure 15(b) illustrates 
an alternative configuration in which the AMCD momentum vector is alined perpendicular 
to the orbit plane for Earth viewing and perpendicular to the ecliptic plane for solar 
viewing. This configuration allows complete spacecraft reorientation about the AMCD 
spin axis with no requirement for direct external torques, orbit regression being 
neglected. The second configuration does require that pointing control torques be gen- 
erated about the spin axis; this was accomplished by using the spin motor with an inherent 
increase in required power. The pointing field of view in both cases is limited at least 
partially by the magnetic bearing gap and thus an additional mechanism is required for 
spacecraft maneuvering. 
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Active magnetic bearing/torquer 



(a) Spin axis pointing. 



(b) Spin axis perpendicular to pointing axis. 


Figure 15.- Active AMCD-spacecraft system. 


AMCD-Spacecraft Maneuvering 

Large angle spacecraft reorientation or maneuvering can be accomplished in a 
variety of ways by using one or more AMCDs; however, only the two ways which appear 
to be most practical for large radius AMCDs are discussed. In figure 16(a), a single - 
gimbal, single AMCD for arbitrary orientation of the pointing axis of a spacecraft is 
shown. By making an analogy with Earth -based telescopes, the azimuth angle a is 
generated by AMCD spin -axis control, and the elevation angle /3 is generated by pro- 
ducing torques against the AMCD momentum on the spacecraft by using a gimbal 
torquer(s). By limiting the spacecraft momentum during the maneuver to a small frac- 
tion of the AMCD momentum, one can see that the AMCD would remain nearly fixed, and 
the spacecraft would move in a fashion similar to that of an Earth-based telescope. When 
the spacecraft pointing axis is appropriately alined with the new target, the gimbal torquer 
would be locked in place and fine pointing would be accomplished by using the active tech- 
nique discussed previously. 
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(b) Counterrotating AMCD pair configuration. 

Figure 16.- Two AMCD- spacecraft maneuver configurations. 

In figure 16(b) a second configuration is shown where two counterrotating AMCDs 
are shown mounted directly to a spacecraft. By reorienting the two rims within their 
respective magnetic bearing gaps or by using a limited gimbaling of one AMCD and also 
spinning up or despinning the rims, some momentum with appropriate direction can be 
imparted to the spacecraft. This momentum represents a maneuvering spacecraft with 
termination of the maneuver begun at a desired orientation. If all initial (prior to the 
maneuver) nonzero momentum can be accommodated by the two AMCDs at the new orien- 
tation, all spacecraft rates can be nulled. Thus, arbitrary maneuvers can be performed 
with this system if given sufficient magnetic bearing gap or gimbal freedom and spin 
speed range. Once the spacecraft has reached a new orientation, one or both of the 
AMCDs can be used in an active or passive mode for spacecraft pointing or stabilization. 
This configuration can also be used for combined power storage and control by maintain- 
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ing a net momentum difference between and H 2 . The net momentum vector can 
stabilize the spacecraft in a passive manner as discussed previously, and power can be 
put into and taken from the system by spinning the two rims up or down while maintain- 
ing a fixed difference in rotational speed. 

Additional Applications 

The AMCD can be used as a vernier experiment pointing system by mounting the 
experiment to the rim and the magnetic bearing to the spacecraft (or vice versa) and 
providing vernier pointing torques with the magnetic bearings by using fine pointing 
information from the experiment package. The rim does not spin at high speed but the 
spin motor is required for roll orientation and control. The magnetic bearings produce 
torques against the spacecraft inertia and stabilization system rather than against the 
rim momentum. 

Although this report has emphasized the large radius applications of the AMCD, 
the device may well be appropriate as the spin assembly for any conventional momentum 
storage device, such as a reaction wheel, control moment gyroscope, or the gyroscopic 
sensing systems found in aerospace vehicles. 

CONCLUDING REMARKS 

A device with potential broad application to spacecraft control and other momentum 
storage applications has been defined and considerations of its configuration and use dis- 
cussed. Examination of the device as a momentum storage unit and system indicates the 
following advantages: 

1. The rotating element approaches a thin rim which is the optimum shape for a 
given stress -limited material when maximizing momentum for a given mass of material 
and for a given maximum radius. 

2. The thin rim allows a unidirectional filament layup of composites; thus, the 
maximum usage of these high strength -weight materials is allowed. 

3. The configuration allows, where possible and desirable, the use of a large rim 
diameter with the inherent additional increase in momentum -mass ratio. 

4. The noncontacting magnetic bearings and drive motor eliminate all mechanical 
friction and wear and should yield a device reliability equal to that of the solid-state 
electronic circuits. 

5. The isolation of the rotating rim from the spacecraft affords an effective control 
over transmittal of rim vibration to the spacecraft when active magnetic bearings with no 
permanent magnetism are used. This is true because any forces transmitted to the 
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spacecraft would be the result of sensed rim motion and subsequent flux coil energizing, 
and would thus allow appropriate notch filtering in the electronics at spin and vibrational 
frequencies. 

6. The magnetic bearings also provide the capability for directly producing torques 
on the spacecraft with no mechanical or electrical breakout torques involved. 

7. For the "passive” mode of spacecraft control (when compared with single-spin, 
dual- spin, or gyrostat control), much improved precessional damping can be shown 
theoretically. 

8. For the "passive” mode of spacecraft control, smaller attitude errors caused 
by environmental torques will result from the higher momentum allowed by the AMCD 
for a given momentum storage weight. 

9. For the "active" mode of spacecraft control, extreme precision fine pointing is 
projected since extremely low spacecraft control torques can be easily generated with 
the magnetic bearings used as spacecraft torquers driven by spacecraft attitude sensors. 

Finally, it is noted that the AMCD can be used as the spin assembly for any of the 
conventional momentum storage control devices (such as reaction wheels or control 
moment gyroscopes) with the unit characteristics of the AMCD offering the potential for 
significant improvements. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., December 19, 1974. 
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APPENDIX A 


FLYWHEEL MOMENTUM DENSITY 


By considering the constraint on flywheel speed to be determined by allowable 
working stress, the specific energy of a flywheel can be determined from the relation- 
ship (ref. 2) 



(Al) 


where Tj^ is the stored kinetic energy; m, the mass; K s , a dimensionless flywheel 
shape factor; a, the design working stress; and p, the material density. 

This relationship (eq. (Al)) can be modified to yield the specific momentum of a 
flywheel by noting that 



H - mk 2 w 


(A2) 


where H is the angular momentum; w, the angular velocity; and k, the radius of gyra- 
tion of the flywheel. 

Substitution of equations (A2) into equation (Al) yields 


§ = k(2K s 2 ) V2 (A3) 

Shape factors for homogeneous materials are given in reference 2 for nine different 
flywheel geometries and vary from 0.305 for a flat pierced disk to a theoretical limit of 
1.000 for a constant -stress disk with an outer diameter approaching infinity. The shape 
factor for a thin rim is 0.500. A factor of 0.85 represents a practical limit. 
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APPENDIX B 


AMCD STIFFNESS EFFECTS 


The equation of motion for AMCD rim vibrations in the spin plane is assumed to be 
that of a rod under tension with rim thickness negligible when compared with rim radius. 
The equation, from reference 8, is 


8S “ 9s^ 8t z 


(Bl) 


where a is the tangential stress; A, the cross-sectional area; Y, Young's modulus; 

I, the area moment of inertia; p, the density of the rim; 77 , the amplitude of in-plane 
vibration; and s, the distance along the rim. The frequency of the fundamental mode 
can be found from the following solution to equation (Bl): 

V - V 0 sin (2 S) sin 0 ) f t (B2) 

where r} 0 is an arbitrary constant. Also it is noted that cr = pr2a;2 f or a thin rim, 
where co is the spin frequency. 

Substituting equation (B2) into equation (Bl) and simplifying yields the relationship 

w f = \J (2w) 2 + (B3) 

where 




(B4) 


From equation (B3) it can be seen that the vibration frequency is always at 
least twice as great as the spin frequency a,’. This formulation neglects rim inexten- 
sibility for the sake of simplicity. 
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APPENDIX C 


DERIVATION OF "PASSIVE” AMCD -SPACECRAFT EQUATIONS OF MOTION 


The derivation begins by considering Euler’s equations of motion in matrix nota- 
tion for an arbitrary rigid body written with respect to a set of axes fixed to the body, 

T b = ^b + w b H b ( C1 ) 

where Tjj are the external torques; Ifo, the inertia matrix; the body rates of rota- 
tion written as a cross product matrix; and H^, the angular momentum. This body 
represents the rotating rim of the AMCD. 

Since the magnetic bearings which produce the torques T^ are fixed to the space- 
craft, it is required to transform equation (Cl) to a second arbitrarily oriented axis sys- 
tem. This transformation E a b is defined by the equation 

v a - E ab v b (C2) 


where v a and are arbitrary vectors. 

The transformation of equation (Cl) begins by multiplying both sides by E a b and 
by substituting, as indicated, the unity matrix taken as E"^E ab . This procedure yields 


E ab T b - E ab I b E ab E ab cI, b + E ab w b E ab E ab H b 


(C3) 


Equation (C3) can be simplified by defining 
~ E ab T b 


Iq = E q i,Ii,E 


-1 


*a “ “abm^ab 
H a = E ab^b 




= E ab aj b E ab 


J 


(matrix) 

w a = E a b w b (vector) 

By differentiating the last of equations (C4), it is seen that 
^a = E ab*b + E ab w b ] 


E ab^b = "a - E ab w b 


Substituting equations (C4) and (C5) into equation (C3) yields 


- 1 . 


Ta = ^a^a + w a E a - ^a E ab E ab w a 


(C4) 


(C5) 


(C6) 
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APPENDIX C - Continued 


By assuming the a-coordinate system and the b-coordinate system to have the same 
Z-axis, with a rotation of | about the Z-axis of the b-coordinate system relative to the 
a-coordinate system, 


cos £ 


~au 


Iffin £ 

i => 


o 


-sin | 0 

cos | 0 

0 1 


(C7) 


Equation (C6) can now be expanded and, by assuming AMCD rim symmetry 
(Ibx - Iby = ! a)» becomes 

Tax “ la “ax + (laz - Ia)“ay“az + ^a“ay 
Tay " ^a “ay + (^a " ^az)“ax“az " ^a“ax / 

Taz = Iaz“az 


<C8) 


Equations (C8) can now be modified to separate the effect of the large AMCD rim 
spin velocity by introducing the variable o>az which will represent the angular velocity 
of the a-coordinate system about the Z-axis where 

“az = “az + £ (C9) 


Equations (C8) now become 

T ax = r a“ax + (!az - J a)“ay“az + !az£“ay 
Tay = *a“ay + (*a " *az)“az“ax - *az^“ax / 
Taz = ^az“az + ^az£ 


(CIO) 


These equations are the equations of motion of the AMCD rim written in the a-system. 

A second set of Euler's equations for the spacecraft is now introduced, where the 
algebra of the transformation of axes is omitted and where the AMCD and spacecraft 
centers of mass are assumed to be coincident, 

Qx + T sx = ^s^sx + (*sz _ I s)“sy“sz + !sz^“sy 

Gy + T S y = Igd>gy + (ig - Iaz)“sz“sx ~ Isz^“sx > (Cll) 

*^Z + Tgz - Isz^sz + ^SZ^ 

J 

where G is an external disturbance torque, spacecraft symmetry is assumed 
(l S x = Isy = Is)? and “sz = “sz + 0. These equations are written in the s -system. 
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APPENDIX C - Continued 


The rates oj£ z and a >§ z are set equal to zero and the a- and s -coordinate axes 
are assumed to be nearly coincident except for small transverse relative rotations. This 
procedure allows a simplification in the calculation of the interaction torques by intro- 
ducing a transformation matrix from the s -coordinate system to the a-coordinate system 
as follows: 

E as = E ai E is = E 2 (0 a ) E 1 (^ a )[E 2 (0 s ) Ejfog)] - 1 = E 2 (0 a ) E 2 (0 a ) Ef 1 !^) E 2 ~1(0 S ) 

= E 2 (0 a ) Ej(tf> a ) Ej(-0 g ) E 2 (-0 g ) = E 2 (0 a ) E^ - <P S ) E 2 (-0 s ) (C12) 


where the subscript i is introduced to represent an inertial reference and the notation 
Ej( ), j = 1 and 2, refers to a transformation from one coordinate system to another 
which has been rotated through an angle denoted by ( ) about an axis j. The Euler 
angles chosen to represent the positions of the a- and s -coordinate axes with respect to 
an inertial axis set are <p and 0 with a 1,2 rotation sequence (x,y) selected. 

By expanding the matrices in equation (C12) for the case when <£ a - $ s , 0 a , 
and 6 S are small, the following equation is found to hold: 

Eas = E 2(0a - ®s) E]/0 a - 0s) = Ej(0 a - 0 S ) E2(0 a - $s) (C13) 


Because of the small angles (<£ a - <p s ) and (0 a - 0 S ), Euler’s theorem on rotations can be 
used to show that 


X = 


(0a * 0 s f + ( 0 a - 



1/2 


sin i 



6 


s 


cos i = 


0a “ 0 


s 




di _ (^a ~ ^s) cos \ ~ (0a ~ 0s) s * n * 
dt ” X 


where the second set of Euler angles X and i is shown in figure 17. 

The torques on the AMCD rim resulting from "passive” or centering-only magnetic 
bearing forces Fj_, F 2 , and F3 (see fig. 17) can be expressed as a function of the rota- 
tion X and the angle y and the linear magnetic bearing gains Kp and Kj as follows. 
The bearing forces acting on the rim are defined by 

Fj = (KgX + K£x)r sin y + KfrXry cos y ^ 

F 2 = -(K^X + K^x)r sin (60° + y) - K jXr y cos (60° + y) \ (C14) 

F 3 = (K^X + K^x)r sin (60° - y) - K^Xry cos (60° - y) J 
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APPENDIX C - Continued 



Figure 17.- Magnetic bearing forces and 
resultant torque. 


where y is shown in figure 17. The sum of Fj, F 2 ,and F 3 is zero. The forces 
are assumed to act perpendicular to the rim. 

The torque is the resultant of the individual torques produced by the bearing forces 


T = [-Fj sin y + F 2 sin (60° + y) - F 3 sin (60° - y)]r 
Tp = [-Fj cos y + F 2 cos (60° + y) + F 3 cos (60° - y)j 



(Cl 5) 


Substituting equations (C14) into equations (C15) and simplifying yields 


T - -K X X - K X X 
T p = -K^Xy 


(C16) 


where 

K X = 1.5r2K,') 

\ (C17) 

= 1.5r 2 K£ j 

Thus, for the case of 120° spacing, the net bearing torque is not dependent on the position 
of the magnetic bearing segments relative to the rim (the torque is independent of y) and 
therefore the three linear bearings may be treated as two rotational bearings. The bear- 
ing torques can now be written as 
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APPENDIX C - Concluded 


T ax - -K x(0a - 0 S ) " K \(0a “ 0s) - K iH e a. " ®s ) 
T ay - -K x (0 a - 9s) “ Kx(®a - ®s) + K X«(0a “ 0s) 


(C18) 


The spin motor torque can be used to control either spacecraft Z-axis attitude or 
attitude' rate as well as to counteract rim drag torque (hysteresis and eddy current 
losses). In this derivation the component of motor torque which is greater than the 
drag torque is of interest and 


Taz - T c - Tfj 

where the subscripts c and d refer to control and drag. 
Note that 


(C19) 


T s - -E sa T a (C20) 

Thus, equations (CIO) to (C12) and (C17) to (C20) represent the nonlinear AMCD- 
spacecraft '’passive” case equations of motion. 


These equations can be linearized by assuming small angular motions except about 
the spin axis, and for the case where T c = T^ and G = 0, they become 


Ia0a - -K\(0 a - 0 s) _ E x(0 a * 0 s) “ E az®a “ K x 6(0 a - 8 s ) 
la® a = _E \(®a ■ ®s) - E \(®a. _ ®s) + E az0a + ^X^0a “ 0s) 

> 

Is0s - K\(0a - 0s) + E \(0a ■ 0s) - H SZ 0 S + K^8($a - 8 s ) 

^s® s = E \(®a _ ®s) + E x(®a - ®s) + ^sz0s ” ^X^(0a 0s) 

J 

where 


(C21) 


^sz _ 


^az “ *az£ 


The linearized equations for the Z-axis represent simple single-axis reaction-wheel 
dynamics and are not considered. 
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